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Abstract. The CELESTE experiment uses the heliostats of an old solar farm in 
the French Pyrenees to detect gamma ray air showers by the atmospheric Cerenkov 
technique. Observations of the TeV blazar Markarian 421 have been made with the fully 
instrumented CELESTE experiment since December 1999. The detection of gamma 
ray emission from this source at energies greater than 50 GeV is presented here. A 
comparison is made with the light curve from the CAT experiment at the same site 
which shows correlation between the observed gamma ray fluxes and the detection of 
short duration flaring episodes. 



INTRODUCTION 

CELESTE has been in stable operation with a 40 heliostat array since October 
1999. Our main targets during this initial period have been the three most reliably 
detected sources in the TeV sky; the Crab nebula, Markarian 421 and Markarian 
501. Our detection of the Crab nebula is described elsewhere in these proceedings 
[1]. This report is concerned with the observations of Mkn 421. 

Mkn 421 is the closest X-ray selected BL Lac object (z=0.030) and was the first 
extragalactic object detected at TeV energies [2]. The emission from this source is 
extremely variable over timescales as short as 15 minutes [3] and multiwavelength 
observations have shown that the TeV variability is correlated with the X-ray emis- 
sion [4]. Observations in the 50 GeV region, not possible until now, will help to 
improve our understanding of the emission processes at work. 

The CELESTE experiment is described in detail in the experimental proposal 
[5] and in [6,7]. Briefly, CELESTE uses forty, 54 m 2 heliostats of a former solar 
electrical plant at the Themis site in the French Pyrenees (N. 42.50°, E. 1.97°, 



altitude 1650 m). The light from the seperate heliostats is resolved using secondary 
optics at the top of a 100 m tall tower. At the secondary mirror focus is one PMT 
for each heliostat. A solid Winston cone at the face of each PMT defines the field 
of view of each tube to be 10 mrad. 

The trigger is constructed by summing eight PMT signals in five groups and 
sending each sum to a discriminator with a threshold set to 4.5 photo-electrons 
(PE) per heliostat. A logic coincidence of at least three of the five groups triggers 
the experiment. Each PMT signal is sent to a 1 GHz flash ADC circuit which 
digitizes the signal. When a trigger occurs, digitization stops and a window of 
100 ns centered at the nominal Cerenkov pulse arrival time is read out. 

ANALYSIS 

CELESTE, with its large mirror area per PMT, is very susceptible to systematic 
effects in the data caused by differences in the night sky background light between 
the ON source observation and OFF source control data. In the case of Mkn 421, a 
star of magnitude 6.1 located at the source position causes a difference of ~15% in 
the PMT anode currents between the two fields. To compensate, we have adapted 
the software padding method proposed by Cawley [8] for use with our FADC data 
[9]. We also apply a software trigger by using the FADC data to reconstruct the 
five trigger group sums and then demanding: >4 groups >5.0 PE per heliostat. This 
increases the energy threshold, but provides us with comparable background data 
in the ON and OFF source fields. We also require > 10 Cerenkov peaks with an 
amplitude > 25 digital counts to ensure the shower parameters are measurable. 

The rest of the analysis is concerned with improving the gamma ray signal to 
hadron background ratio using two parameters. The first, <J grp , is a measurement of 
the homogeneity of the light distribution between the trigger groups and should be 
the most efficient background rejection cut according to our simulations, rejecting 
80% of the remaining hadrons for only 30% of the gammas. The second parameter, 
9, is the shower axis angle relative to the pointing direction reconstructed using 
two points; the barycentre of the light distribution at ground level and the position 
of the shower core relative to the pointing direction at the altitude of maximum 
emission deduced from a sphere fit to the Cerenkov pulse arrival times at the 
heliostats. 

OBSERVATIONS AND RESULTS 

Mkn 421 was observed between November 1999 and March 2000. After rejecting 
those observations with bad weather or equipment problems, the dataset consists 
of 80 ON/OFF pairs with an ON source exposure of 22.5 hours. 

The CAT imaging telescope [10], situated on the same site as CELESTE, made 
observations of Mkn 421 during the same period. In December 1999, no evidence for 
a signal was reported by CAT. The left-hand plot of Fig. 1 shows the distribution 
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FIGURE 1. The difference in the distribution of the shower axis angle 9 between ON and OFF 
source observations of Mkn 421. Left: For the December 1999 quiet period. Right: For the 
January/February 2000 active period. 



of the difference in the shower orientation parameter, 6, between the ON and OFF 
source observations as measured by CELESTE during December. No evidence for 
a signal at low 9 is seen (for 1.5 hours of ON source exposure). Table 1 shows the 
number of events remaining at each stage of the analysis. 

In January/February 2000 CAT reported that Mkn 421 was in an active state. In 
particular, two bright flares were observed, with emission levels >2 Crab. The right- 
hand plot of Fig. 1 shows the difference distribution in 8 for CELESTE observations 
during these flares (5.2 hours of ON source exposure). A strong signal is apparent 
at low values of 6>, as expected for gamma rays. The significance of the excess is 
8.1 a. Table 2 shows the improvement in the signal/background ratio at each stage 
of the analysis after the initial cleaning cuts. The significance before cuts is not 
very meaningful as the excess, while undoubtably largely due to gamma rays, is 
contaminated by systematic effects due to sky noise differences. 

The whole Mkn 421 dataset has been analysed and Fig. 2 shows the lightcurve 
as observed by both CAT and CELESTE on a night by night basis. Emission 
levels up to ~ 3 Crab (~ 7 7 min -1 ) were observed, detectable by CELESTE in a 
single night. Also shown in the figure is a plot of the gamma ray flux measured by 
CELESTE against that measured by CAT showing reasonably good correlation. 



TABLE 1. The number of events, excess, significance and signal to back- 
ground ratio (S/B) after each cut for the December quiet period. 





ON 


OFF 


Excess 


Significance 


S/B 


Raw 


104443 


10 3991 


452 


1.0 


0.4% 


Software Trigger 


55 294 


55 236 


58 


0.2 


0.1% 


iV poaks > 10 


51943 


51924 


19 


0.1 


0.0% 


(Tgrp < 0.25 


12154 


12186 


-32 


-0.2 


-0.3% 


9 < 7 mrad 


6 648 


6 679 


-31 


-0.3 


-0.5% 



TABLE 2. The same as Table 1 for the January/February active period. 





ON 


OFF 


Excess 


Significance 


S/B 


7/min 


Raw 


369 424 


362 248 


7176 


8.4 


2.0% 




Software Trigger 


177 321 


174 273 


3 048 


5.1 


1.7% 


9.8 


Speaks > 10 


167 385 


164 345 


3 040 


5.3 


1.8% 


9.8 


tTgrp < 0.25 


29 049 


27 536 


1513 


6.4 


5.5% 


4.9 


9 < 7 mrad 


17169 


15 709 


1460 


8.1 


9.3% 


4.7 





X /ndf = 15.7/13 



FIGURE 2. The CAT/CELESTE night by night light curve for Mkn 421 and the correlation 
between the experiments for measurements taken on the same nights. 



CONCLUSIONS 

The TeV blazar Mkn 421 has been detected by CELESTE for the first time at 
an energy of ~ 50 GeV. Evidence for night by night variability has been observed, 
and the variablity measurements by CELESTE are correlated with those of CAT. 
Importantly, a non-detection of the source in a quiescent state, followed by a very 
significant detection in a flaring state demonstrates that the CELESTE analysis is 
working, and that we have been able to correct for systematic effects due to night 
sky light differences between the ON and OFF source fields. 
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